Optimal control of entanglement via quantum feedback 
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It has recently been shown that finding the optimal measurement on the environment for sta- 
tionary Linear Quadratic Gaussian control problems is a semi-definite program. We apply this 
technique to the control of the EPR-correlations between two bosonic modes interacting via a para- 
metric Hamiltonian at steady state. The optimal measurement turns out to be nonlocal homodyne 
measurement — the outputs of the two modes must be combined before measurement. We also find 
the optimal local measurement and control technique. This gives the same degree of entanglement 
but a higher degree of purity than the local technique previously considered [S. Mancini, Phys. Rev. 
A 73, 010304(R) (2006)]. 

PACS numbers: 03.67.Mn, 02.30.Yy, 42.50.Lc, 42.50.Dv 



> 

o 
o 
o 

\o ■ 

o ; 

^ ■ 

a,: 

^— > ■ 



X 



I. INTRODUCTION 

Quantum feedback control is a well-established theo- 
retical technique for stabilizing an open quantum sys- 
tem in a state with certain desired properties 0, 0, 
The basic idea is to use the information that leaks from 
the system into a bath to undo the undesirable effects 
of coupling to this, or other, baths. Notable examples 
include protecting a "Schrodingcr cat" superposition 14, 
correcting errors in encoded quantum informational, |6| , 
maintaining a two-level atom in an arbitrary state pjhde- 
terministically producing entangled states of s pins M 
(which has been experimentally demonstrated and 
cooling various systems to (close to) their ground states 

IH El EH. 

Recently, one of us started to consider the applica- 
tion of quantum feedback to control entanglement. Pre- 
liminary studies have been carried out for two interact- 
ing qubits H and two interacting bosonic modes 0], 
damped to independent baths. Physically, two damped 
and interacting bosonic modes could be realized by opti- 
cal cavity modes coupled by a %' 2 ' nonlinearity. Such a 
nonlincarity induces only a finite amount of entanglement 
between the modes in steady-state. By contrast, it was 
shown in Ref. H] that performing homodyne detection 
on the two outputs, and using these currents to modulate 
the (linear) driving of the two modes could, under ideal 
conditions, increase this entanglement without limit. 

The quantum control problem of Ref. 1 161 has, like 
many which have been considered 0,0jI11jHE1> an anal- 
ogy in the class of classical LQG problems H- That is, 
systems with Linear dynamics and a Linear map from 
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inputs to outputs, an aim that can be expressed in terms 
of minimizing a Quadratic cost function, and Gaussian 
noise in the dynamics and the outputs. The quantum 
LQG problem has recently been analysed in detail by 
one of us and Doherty |l8j . In particular, it was shown 
in Ref. |l8j that in the quantum case, there is an extra 
level of optimization that naturally arises: choosing the 
optimal unravelling (way to extract information from the 
bath) given a fixed system-bath coupling. 

In this paper we reconsider the problem of Ref. H 
from the perspective of Ref. ■ We formulate the prob- 
lem as an LQG control problem, and find the optimal 
unravel ling . This is different from the unravelling used 
in Ref. jlfj (it requires interfering the two output beams 
prior to homodyne detection) and leads to greater entan- 
glement for any strength of the nonlinearity. This shows 
the usefulness of the general techniques of Ref. [T^j . 

The paper is organized as follows. In Sec. [H] we sum- 
marize the general theory of quantum LQG control prob- 
lems from Ref. |18| as needed for the current problem. 
In Sec. IIIII the problem of maximizing the steady-state 
entanglement between two bosonic modes interacting via 
parametric Hamiltonian is addressed and the optimal un- 
ravelling found. As stated above, this involves interfer- 
ing the output beams from the two modes prior to detec- 
tion, which could be regarded as a nonlocal measurement. 
In Sec. IIVI we consider the constraint of local measure- 
ments (that is, independent measurements on the two 
outputs). We consider a variety of measurement and 
feedback schemes, including homodyne and heterodyne 
detection, and find that two schemes, that considered in 
Ref. |l6j| , and another (more symmetric) scheme, are the 
best. However, when it comes to the purity of the sta- 
tionary entangled state, considered in Sec[V] the more 
symmetric scheme is superior. Finally, Sec I VII is for con- 
clusions. 
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II. QUANTUM LQG CONTROL PROBLEMS 

A. Continuous Markovian Unravellings of Open 
Systems 

We consider open quantum systems whose average 
(that is, unconditional) evolution can be described by an 
autonomous differential equation for the state matrix p. 
The most general such equation is the Lindblad master 
equation 



P = -i[H, p] + T>[c]p = C p 



(1) 



Here H = W is the system Hamiltonian (we use h = 1 
throughout the paper), while c is a vector of operators 
c = (ci, • • • , cl) t that need not be Hermitian. The action 
of D[c] on an arbitrary operator p is defined by 



V[c]p^J2 



i=i 



cipc\ - - (cjcip 



pcjci 



(2) 



Master equations of this form can typically be derived if 
the system is coupled weakly to an environment that is 
large (i.e. with dense energy levels). Under these condi- 
tions, it is possible to measure the environment continu- 
ally on a time scale much shorter than any system time 
of interest. This monitoring yields information about the 
system, producing a stochastic conditional system state 
p c that on average reproduces the unconditional state p. 
That is, the master equation is unravelled into stochastic 
quantum trajectories ppl ]. with different measurements 
on the environment leading to different unravellings. 

For the purposes of this paper we can restrict to un- 
ravelings that yield an evolution for p c that is continuous 
and Markovian. In that case, it must be of the form [2lJ 

dp c = Capcdt + dz\t)A c cp c + p c A c ctdz(t). (3) 

Note that here the f indicates transpose (T) of the vector 
and Hermitian adjoint of its components. We are also us- 
ing the notation A c 6 = 6 — (6) , where (6) = Tr[p c o]. Fi- 
nally, we have introduced a vector dz = (dz u • ■ ■ ,dzL) T 
of infinitesimal complex Wiener increments j22|. It sat- 
isfies E[dz] = 0, where E denotes expectation value, and 
for efficient detection has the correlations 21] 



dzdz^ = Idt, dzdz T = Tdt. 



(4) 



Here T is a symmetric complex matrix, which is con- 
strained only by the condition U > where U is the 
unravelling matrix 



I + Rc [T] 
Im [T] 



Im [T] 
I - Re [XI 



(5) 



Equation describes a quantum diffusion process for 
the conditional state p c . Equations of this form were 
first written down by Belavk in Jj| , and were derived in- 
dependently by Carmichael [20( to describe homodync 



detection in quantum optics. The measurement results 
upon which the evolution of p c is conditioned is a vector 
of complex functions 



3 T dt = (c T +c t T) c dt + rfz 1 



(6) 



Following the terminology from quantum optics, we will 
call J a current. 



B. Linear Systems 

We now specialize to systems of N degrees of freedoms, 
with the nth described by the canonically conjugate pair 
obeying the commutation relations [q n ,Pn] = i- Defining 
a vector of operators 



Pi, -,qN,PN) 



(7) 



we can write 



[%m — ^nmi ($) 

where S is the (2N) x (2N) symplcctic matrix 



A' 



E -et-°ii 

n=l v 



-s- 1 . 



(9) 



For a system with such a phase-space structure we can 
define a Gaussian state as one with a Gaussian Wigner 
function . We write the mean vector as (x) and its 
covariance matrix as V: 



((Ax n Ax m ) + (Ax m Ax n )) j2. 



(10) 



For these to define a quantum states, the necessary and 
sufficient condition is that I2J1 



U + iS/2 > 0. 



(11) 



To obtain linear dynamics for our system in phase- 
space we require that H be quadratic, and c linear, in 



H = (l/2)x T Gx-x T EBu(t), c=Cx, 



(12) 



where G is real and symmetric and B is real. The second 
term in H is linear in x to ensure a linear map between 
the time-dependent classical input u(t) to the system and 
the output current J(i). For such a system, the uncon- 
ditional master equation has a Gaussian state as its 
solution, with the following moment equations 



d{k)/dt = A(x) + Bu(t), 
dV/dt = AV + VA T + D. 



(13) 
(14) 



Here A = S(G + lm[&C}) and D = SRc[C' t C , ]S T . 

For conditional evolution of linear quantum systems it 
is convenient to recast the complex current J of Eq. (jSJ) 
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as a real current with uncorrelated noises: as opposed to 
the complex current J with (in general) correlated noises: 



Here C = 2([/) 1 / 2 C, where 

C T = (Re[C T ],Im[C* T ] 



(15) 



(16) 



while dw is a vector of real Wiener increments satisfying 
dwdw T = Idt. For linear systems this conditional state 
p c from Eq. J2J) is Gaussian, with the conditional moment 
equations [18j 

d(x) c = L4(x) c + Bu(t)} dt + (V C C T + T T ) dw (17) 
V c = AV C + V C A T + D- {V C C T + T T )(CV C + T). (18) 

Here T = -(Uy/ 2 SCE T , where 

/ 



S = 



-I 



(19) 



Note that the equation for V c is deterministic. In many 
situations (including those considered later in this pa- 
per), Eq. l|18|l has a unique steady-state solution [lSl ] 
which moreover is a so-called stabilizing solution |25j . We 
will notate such a solution V^ s as Wjj to emphasize that 
it depends upon the unraveling U: 



= flWu + Wufl T - WuC T CWu + EE 1 



(20) 



Here f2 = S[G + C T S(2U - I)C] is Hamiltonian drift, 
while E = SC T /2 manifests the measurement back- 
action noise resulting from having y oc Cx. The set of 
Wjjs, for all possible unravellings U (including ineffecient 
monitoring) , is the set of real symmetric matrices satis- 
fying the two linear matrix inequalities Wjj + iS/2 > 
and 



D + AW V + WuA T > 0. 



(21) 



This is shown in Ref. [l8j- Moreover, given a Wu that 
satisfies these inequalities, a (not necessarily unique) un- 
ravelling U that will generate it can be found by solving 



R T UR = D + AWu + WuA T , 



(22) 



where R = 2CWu/ + SCE. 



C. Optimal Quantum Control 

In feedback control, u(i) depends on the history of the 
measurement record y(s) for s < t. The typical aim of 
control over some interval \tg, tA is to minimize the ex- 
pected value of a cost function |l7| . the integral of the 
sum of positive functions of x(t) and u(t) for to < t < t\. 
We are interested in the special case of Linear-Quadratic- 
Gaussian (LQG) control [l7T |: a linear system with a 



quadratic cost function and having Gaussian noise. For 
an LQG control problem it can be shown that the optimal 
u is linear in the phase-space mean: 



u(t) = -K{t) (x) c (i), 



(23) 



where the matrix K (t) can be determined from A, B, 
and the cost functions. It is independent of D, C, and T. 

In this paper we arc concerned only with the properties 
of the system at steady-state, so our aim is to minimize 
to = E[/i] in the limit t\ — * oo, where 



/> = <x T Px) c , 
with P > 0. Note that in steady state 



(24) 



E ss [(x T Px)J = tr[WuP] +E ss [(x) c T P(x)J. (25) 

Assuming (as is the case in our system) a stabilizing solu- 
tion Wu plus control over all relevant degrees of freedoms 
of the system (as will be the case if B is invcrtible), the 
control can always be chosen to set (x) c — > 0, so that 

TO opt =tt[PWu]. (26) 

For such systems it turns out 01 that the same result 
(that is, (x) — ► 0) can always be achieved with Marko- 
vian feedback as introduced by Wiseman and Milburn 
p^ . This is a much simpler form of feedback; for a gen- 
eral linear system it means 



u(t) = F(t)y(t) 



(27) 



If F is time- independent, the average evolution of the 
system is described simply by modifying the drift and 
diffusion matrices to 



A' = A + BFC, 
D' 



(28) 



D + BFF T B 



BFT + T T F T B T . (29) 



With B invertible it can be shown that the optimal choice 
(which makes (x) c -> 0) is BF = -W V C T - T T . 



III. CONTROLLING ENTANGLEMENT 

We now specialize to the system examined in Ref. ^(| : 
a non-degenerate parametric oscillator p7T | where two 
damped bosonic modes c\ and C2 interact through a x' 2 ' 
optical nonlinearity. Treating the pump mode classically, 
this results in a quadratic Hamiltonian for the two modes. 
The master equation is 



p = -i 



H,p 



where 



H = ix [c\c\ - cic 2 



V[ Cl \p + V[c2] P , (30) 



x(<ZiP2 +92Pi) • (31) 



Here \ is the coupling constant, proportional to the x' 2 ' 
coefficient and the amplitude of the pump. We have also 
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defined quadratures for the two modes via Cj — (q~j + 

This system fits within the general model described 
above, with N = L = 2. Defining 



we have 



G 



= (ft, Pi, 92, P2) 



/o x 
0x0 
0x00 
\x 



and 



1 j 00 
001 i 



(32) 



(33) 



(34) 



From the above theory we obtain 



V2 



10 
10 
10 
1 



(35) 



.4 







(-\ 
-± 



X 

V 



D = - 
2 







x \ 

-x 
-1 





(36) 



/ 1 
10 
10 

Vo 1 



(37) 



For x < 1/2, the system has a stationary state with 
mean zero and covariance matrix V which can be found 
by setting dV/dt to zero in Eq. Q14[l. The result can be 
written in terms of 2 x 2 submatrices as 



V 



7 a 
T 

a 7 



where the matrix elements of 7 and a are 



1 





= Ipp = 


2 


Iqp 


= Ipq = 


0. 


V qq 






Pqp 


= V-pq = 


0. 



1 



1 -4x 2 
X 



I -Ax 1 



(38) 

(39) 

(40) 
(41) 
(42) 



Since the steady state is Gaussian it is completely char- 
acterized by the correlation matrix (|38|l . Then, its degree 
of entanglement can be quantified by means of the loga- 
rithmic negativity [2^| 



Ls -log(2C- 



if C < 1 
otherwise 



(43) 



where 



C- 



(det 7 — det a) — y (det 7 — det a)" — dct V 



(44) 

is the lowest symplectic eigenvalue of the partial trans- 
posed Gaussian state characterized by V. The quantity 
L is represented by the curve d) in Fig^ It is nonzero 
for all x > 0, and is finite even as x ~~ * 1/2- That is, 
the damping channels degrade the system state, prevent- 
ing it from becoming maximally entangled like that of 
Ref. HlJ. 

A Gaussian with a covariance matrix of the form of 
Eq. (|38|) is entangled if and only if the variance of the 
mixed quadratures 



(45) 



is less than the the vacuum fluctuation level of unity |2£ 
l30l | . It is easy to calculate that 



((£^9) + x 2 (tt - 6)Y 



1 



l + 2x 



(46) 



We see that the variances l|46f) go below 1 as soon as 
X > 0, from which we infer entanglement. Since we must 
have x < 1/2, the variances l|46(l are limited from below 
by 1/2. This is even though the variance in Xj(9), for 
either j and for all 9, are unbounded as x -* 1/2. This 
again shows that the stationary state has only a finite 
amount of entanglement, and is not pure. For states of 
this form, the log-negativity is in fact a simple function 
of the above variance: 

L = -log 2 [l/(l + 2 X )]. (47) 



A. Optimal Measurement and Control 

As we have seen, entanglement is manifest in the 
squeezing of the quadratures in Eq. I|46|l. for all 9. Thus, 
as an aim for the feedback, we can choose the minimiza- 
tion of 



d9 

-(( Xl (9)+x 2 (n-9)f) 



in steady state. This evaluates to 

((q 1 -q 2 j 2 )/2 + ((p 1 +p 2 ) 2 )/2 
This is exactly of the form of Eq. (|24(l , with 



(48) 



(49) 




(50) 



From the symmetry of the problem, we can assume 
that the optimal conditional covariance matrix shares the 
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same structure as the unconditional matrix, namely: 



Wu 



/ a 

a -13 

(3 a 

\ -/3 a 



Thus the quantity to be minimized is 

m = Tr [PWu] = 2 (a - #) 



(51) 



(52) 



We have to find the minimum of m constrained by 
Wu + («/2)S > and D + AWc; + W V A T > 0. In terms 
of a and /3 this becomes 



min (a — (3) 



a- -yr + 4/3 2 > ; 



--(a±/3)(lT2 X ) >0. 



(53) 
(54) 

(55) 



Taking a = i^l + 4/? 2 we get m = 2 (y 1 + 4/? 2 - 2/?) , 

which decreases monotonically with (3. But from the con- 
dition (|55|l we obtain 



That is, 



/3<X 



Jf v /l + 4/34±2/3 N )(lT2x)>0 



i-x 

1-2 X 



0<X<- 



(56) 



(57) 



Thus, choosing (3 = X jz^ an d a = 5\A + 4/3 2 , we ob- 
tain the minimum m op t- In this case the logarithmic 
negativity takes a simple analytical form 



L 



opt 



log 2 [m opt ] = - log 2 [l - 2x] • 



(58) 



This is represented by the curve a) of Fig^ Note that 
this is unbounded as x —> 1/2. 

If now we wish to know how to achieve this optimal 
result, we can use Eq. l(2*2^l to get 



U 




(59) 



Since U = U 1 / 2 we can easily derive the matrix C of 
Eq. Ull) 



C 




(60) 



This tells us that the optimal unraveling is the measure 
of q\ — q*2 (first and second rows) and p\ + p2 (third and 
fourth rows). Intuitively this makes sense, as it is the 



L 




b) 



O 
d) 



o.i 



0.2 



0.3 



. 4 



X 



0.5 



FIG. 1: The logarithmic negativity L of the steady-state 
quantum state of the non-degenerate OPO, versus the optical 
nonlinearity strength \, for: a) optimal (nonlocal) measure- 
ment and feedback of Sec. Ill A; b) optimal local measurement 
(homodyne) and feedback of Sec. IV A; c) non-optimal local 
measurement (heterodyne) and feedback of Sec. IV B; and d) 
no feedback. 



S1 




S2 





w 



FIG. 2: Schematic representation of feedback action based 
on nonlocal measurements. SI and S2 are the two interacting 
systems and M is a common measurement box. 



variances of these quantities that we wish to minimize, 
from Eq. I|49|l . Note however that these are nonlocal mea- 
surements in the sense that they involve combinations of 
observables belonging to the different subsystems. That 
is, the output beam from mode 1 must be mixed at a 
beam splitter with the output beam from mode 2, and 
then the two beam-splitter outputs subject to homodyne 
detection. One of these detections can measure the out- 
put quadrature corresponding to qi — q2, and the other 
can measure that corresponding to pi+p2- This situation 
of nonlocal measurements is schematically represented in 
FigH 



IV. LOCAL FEEDBACK ACTION 

Having shown that the optimal control protocol in- 
volves a nonlocal measurement, it is natural to ask how 
much improvement this offers over protocols involving 
only local measurements on the two subsystems. The 
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latter would be easier to implement experimentally, es- 
pecially if the two modes had significantly different fre- 
quencies. In previous work |l6j |. only local measurements 
were considered. 



A. Single quadrature measurements 

We begin by considering a homodyne measurement of 
the output beam of each mode. From Eq. i|46|) we see 
that there are no preferred quadratures to be measured 
provided that their angles sum up to tt. Without loss of 
generality we can assume to measure q\ and q% . In terms 
of the parameters of Sec. Ill Al we require 



T = diag(l,l), 



(61) 



so that Ji oc qi and J2 oc 52- 

The Hamiltonian term — yi T 'EBFy would represents 
the feedback Hamiltonian. Since we measure q\ and qi 1 
it is natural to act on the quadrature to q\ — q\ , in order 
to minimize its variance. That is, we choose the feedback 



J 



to be proportional to the conjugate quadrature: 



A- [Ji(t) 



M*)] x [pi -P2] 

M*)] x \pl +V2] • (62) 



Here A± represents possible feedback strengths. Equa- 
tion (|62|) represents the most general feedback action 
that accounts for the symmetry between the two sub- 
systems. Note that this feedback can be performed lo- 
cally because the Hamiltonian l|62|) contains no products 
of operators for both subsystems. However, in general it 
requires classical communication, so that the controller 
for mode 1 can apply a Hamiltonian proportional to J2, 
and vice versa. Eq. I|62|) is obtained by choosing the feed- 
back driving like 



A + + A_ A+ - A_ 

J_ I 

" V2 A+ - A- A+ + A- 





(63) 



As consequence of feedback action the matrices A and D are modified according to Eq. I|28|) and Eq. 11291) to 





A' = 



( + +A H 


X - A- + A+ 

V 



X - A- + A + 
\ 
-i + A_+A + 
-X 



\ 

~x 



-u 



(64) 




/ (1 - A- - \+y + (A- - X+Y 2(1 - A- - A+)(A_ - A 

D= 2\ 2(1 -A_ -A+)(A_ -A + ) (1-A_ - A+) 2 + (A_ - ^ ^ 2 n I ■ 
V 

The stationary covariance matrix that results from these is of the form of Eq. I|38|) with 

-1 + 4(1 + X )A+ - 2(1 + 2 X )A 2 + + A 2 _(-2 + A x + 8A + ) - 4A_(-1 + X + 4A+ - 2A^) 
lqq 2(1 + 2 X -4A_)(-1 + 2 X + 4A + ) ' W 

Iqp = Ipq = , (67) 
-( — 

A 2 _(l - 4A+) - A 2 + + 4A_A^ + x (-l + 2A_ - 2A 2 _ + 2A+ - 2\ 2 + ) 

(1 + 2 X -4A_)(-1 + 2 X + 4A + ) ' (59) 

Vqp = CTpg = , (70) 

"» = -rr^- (71) 



7 W = o ( ) . (68) 



We have maximized the logarithmic negativity Ij43(l constraints are 
over A_|_ and A_. with the constraints that V be a stable 

solution to Eq. ijTlfr. In the range < X < 1/2, these \ ± < I zp (72) 



4 2 



We summarize the results by distinguishing four limit 
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K 




< 

S2 






I 

\M2j 



FIG. 3: Schematic representation of feedback action based 
on local measurements and no classical communication. 5*1 
and 5*2 are the two interacting systems and Ml, M2 axe local 
measurement boxes. 



K 




' * 

S2 







feedback Hamiltonian coincides with that used in 
Ref . 0| . The optimal value of the feedback param- 
eter is A = \, and gives rise to a great improvement 
in the logarithmic negativity with respect to the no- 
feedback case [it corresponds to curve b) of FigQ]. 
By approaching the instability point x ~ * 1/2 the 
logarithmic negativity increases indefinitely. 

iv) If we set X± = ±A we once again require classical 
communication (see Fig0J). The optimal value of 
parameter is A = — x, and gives rise exactly to the 
same values of the logarithmic negativity as for the 
case iii) [thus corresponding to curve b) of FigQ] 
too]. However, in Sec. it will become clear that 
the case iv) is superior to the case iii) in other 
ways. 

That case iv) gives the best result is not surprising, 
since it gives rise to a feedback Hamiltonian that resemble 
that in Eq. i|31[l . once it is remembered that J\ oc q\ and 
J2 oc <72- Note that although in the cases iii) and iv) the 
entanglement increases without bound as x ~ * 1/2, the 
log-negativity is still below that of the optimal nonlocal 
feedback for all values of x as shown by Fig^ 



FIG. 4: Schematic representation of feedback action based 
on local measurements supplied by classical communication. 
51 and 52 are the two interacting systems and Ml, M2 are 
local measurement boxes. 



cases for which L becomes dependent on a single param- 
eter. 

i) If we set X± — A we have a purely local feedback, with- 

out classical communication. This is because the 
first (respectively second) current is used to con- 
trol the first (respectively second) subsystem (see 
FigyJ. This case does not show any improvement 
with respect to the no-feedback case; that is the op- 
timal value of parameter is A = [it corresponds to 
curve d) of Fig^. This is because, with local mea- 
surements and no communication, the correlations 
between the two subsystems cannot be increased. 

ii) If we set A_ = and A + = A we do require classical 

communication (see FigQJ. However, also this case 
does not show improvement with respect to the no- 
feedback case; the optimal value of parameter is 
A = [it corresponds to curve d) of Fig^. This is 
because the corresponding feedback Hamiltonian is 
not effective acting on the antisqueezed quadrature 
Qi + <72- 

iii) If wc set A_ = A and A + = we again require clas- 

sical communication (see Fig0J. In this case the 



B. Joint quadratures measurements 

Since Eq. 1)49(1 contains both q and p, one might think 
that performing joint quadratures measurements in both 
subsystems would be an effective route to controlling en- 
tanglement. Of course it is not possible to measure both 
quadratures with perfect efficiency, but it is possible to 
measure each quadrature with an efficiency of 1/2. This 
can be achieved by heterodyne measurement, for example 
j^. In terms of the parameters of Sec. Ill Al we require 
T = so that Ji oc ci and J2 oc C2. 

Bearing in mind the results of preceding subsection 
(that is, that scheme iv) performed best) we restrict our 
consideration to feedback that gives rise to Hamiltonian 
resembling the one in Eq. l)3ip. Hence we choose the 
feedback driving as 



BF 



( fi 

-fj, 

/] 

V -n 



(73) 



corresponding to the feedback Hamiltonian 



Hn, = -i 



■ I' 



[ji(t)c 2 - </rw4) + (j 2 (i)cr - jj(t)gt) 

(74) 

Here jj, represents the feedback strength and i] = 1/2 
accounts for the half unit efficiency. 
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As consequence of feedback action, the matrices A and D are modified according to Eq. (|28f) . Eq. i|29[l to 

\ 

-x - fi 





A' 



( -| x + M 
-\ 
X + M o -| 

V o -x-n o 




_ 1 

2 



(75) 



+ 2 M 2 





-2/i 








1 + 2/x 2 





2/' 


-2/i 





1 + 2 A t 2 








2/i 





1 + 2/t 



D 7 



Proceeding as above, the stationary covariance matrix elements resulting are given by 

-l + 4xAi + 2/i 2 



799 2(-l + 4( X + M) 2 ) 

"fqp "fpq i 



7pp ■ 



(76) 



(77) 
(78) 



'ip 



X + 2 X /i 2 + 2^ 3 
-1 + 4( X + A*) 2 

&pq = , 



(79) 
(80) 



We have maximized the logarithmic negativity Ij43(l 
over /I with the constraints that V be a stable solution 
to Eq. In the range < x < 1/2 these are 

1 1 

-7,-* < V<7,~X- ( 81 ) 

We summarize the results hereafter. 



v) The optimal value of the parameter /x is found to be 
li = (-1 - 2x + yf\ + 4x 2 )/2. This gives rise to 
a small improvement in the logarithmic negativity 
with respect to the no-feedback case — it corre- 
sponds to curve c) of Fig^ 

Although this case does improve entanglement, it is not 
as good as the best homodyne scheme iv). This can be 
understood as follows. In controlling a quantum system, 
one has always to reaach a tradeoff between information 
gain and introduced disturbance. Heterodyne detection 
allows us to gain information about both system quadra- 
tures, in contrast to homodyne detection, at expenses of 
introducing more noise via the feedback. In our system, 
it is apparent that a high degree of entanglement can be 
produced by controlling only one pair of quadratures, so 
the noise introduced by heterodyne-based feedback pro- 
duces inferior performance relative to homodyne-based 
feedback. In other contexts (with other Hamiltonians) 
heterodyne-based feedback may outperform homodyne- 
based. 



V. PURITY 

The fact that for optimal nonlocal feedback, and lo- 
cal feedback of cases iii), iv), the entanglement can in- 
creased without bound, means that feedback is able to 
recycle the information lost by the system into environ- 
ment through the amplitude damping. However the EPR 
correlations [3lJ imply not only an arbitrarily entangled, 
but also a pure state. We now check what the purity of 
our stationary state is under the various feedback control 
schemes. 

For a Gaussian state the von Neumann entropy can be 
written as [12 

S = g(C+) + g(C-) , (82) 



where 




(84) 

are the symplectic eigenvalues of the Gaussian state char- 
acterized by V. 

We have numerically evaluated the quantity l)82|l for 
nonlocal and local feedback action and the results are 
shown in Fig[SJ The curve a) corresponds to the worst 
cases i), and ii) also corresponding to the no- feedback 
action. Below is the curve b) corresponding to the case 
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a) 



S 




0.1 0.2 0.3 0.4 0.5 

X 

FIG. 5: The von Neumann entropy S of the steady-state 
quantum state of the non-degenerate OPO, versus the optical 
nonlinearity strength \, for: a) no feedback; b) non-optimal 
local measurement (homodyne) and feedback of Sec. IV A; c) 
non-optimal local measurement (heterodyne) and feedback of 
Sec. IV B; d) optimal (nonlocal) measurement and feedback 
of Sec. Ill A and optimal local measurement (homodyne) and 
feedback of Sec. IV A. 



trol steady state EPR-corrclations for two bosonic modes 
interacting via parametric Harniltonian oc \- Both these 
actions allow one to produce arbitrary amounts of entan- 
glement as x - * 1/2, although more in the former case. 
Moreover, they both do this while producing a pure state 
- that is they permit us recover the coherence of our 
open quantum system. (Incidentally the possibility of 
coherence recovery by means of feedback was forecast in 
Ref. |34j for finite dimensional systems by an information 
theoretic approach.) 

Our local feedback action requires only classical com- 
munication and Gaussian operations (linear displace- 
ments). This may appear to contradict the impossibility 
to enhance (distill) entanglement by means of Gaussian 
LOCC stated in Refs.j35|. The key point is that, in con- 
trast with Refs.j3f|, here the LOCC operations contin- 
uously happen while the entangling interaction is "on" . 
Thus, the presented approach may shed some light on 
the subject of entanglement distillation. 



iii) and showing that the state remains not pure. In 
this case we have the mixedness, as well as the amount 
of entanglement, increasing by increasing x- The fact 
that they both increase indefinitely may sound strange. 
However, the limit x — > 1/2 has to be taken with care, 
and the above results are justified by the fact that it 
allows infinite energy to come into the state. The curve 
c) corresponds to the case v). 

Finally, the curve d) corresponds to optimal nonlocal 
feedback and case iv), thus showing that in such cases 
the entropy is always zero and the purity of the state is 
restored by the feedback action. These results (FigJSJl to- 
gether with those of logarithmic negativity (Fig^J clearly 
show the optimality of the feedback scheme iv) among 
the local schemes, and the global optimality of the non- 
local scheme. 



While we have used a semidefinite program to find 
the optimal measurement and feedback action for gen- 
eral LQG systems, to find the optimum local scheme we 
used simple optimization informed by the symmetries of 
the system. The question of defining an efficient pro- 
gram to find the optimal local measurement for feedback 
control of general LQG systems remains open. 
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